Shell filling in non-linear magneto-tunneling spectroscopy of vertical quantum dots 
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We report on non-linear magneto-tunneling experiments carried out in single GaAs vertical quantum 
dots. We show that conduction at low voltage bias can be a spectroscopic tool for both the ground 
state and first excited states of few electrons on the dot. Increasing voltage a large resonant peak is 
observed and attributed to tunneling across the quasi-continuum of higher excited states. 
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Vertical quantum dots (QDs) obtained on double 
barrier resonant tunneling heterostructures better than 
those realized on lateral 2DEG systems allow to investi- 
gate the few electrons limit by starting from an empty 
dot and progressively charging it. However linear trans- 
port experiments require the addition of a lateral gate 
around the dots: the fabrication of such devices is a real 
technical difficulty and only few realizations have been 
reported in the literature till now |J . 

We report here on new results in electron spectroscopy 
of more conventional two-contacted vertical QDs without 
a lateral gate. These structures are relatively easier to 
fabricate , but the theoretical analysis is difficult because 
transport is non-linear. 

Coulomb staircase has been observed in these struc- 
tures up to now only in presence of a strong barrier- 
thickness asymmetry |^||: when the electrons enter 
the QD from the thin barrier side, charge accumula- 
tion with subsequent electron-electron interaction is en- 
hanced. Viceversa, one-electron states have been ob- 
served when electrons enter the QD from the thick barrier 
side 0. We have used devices with symmetric barriers 
and, as discussed in the following, the conductance pat- 
tern at low bias can be attributed both to the ground 
states and to the excited states for N-electrons (N < 6). 
At larger voltage bias a crossover to a quasi continuum of 
excited states in the dot produces a huge resonant peak 
as in diodes which are not laterally confined. 

The double barrier structures were grown by molecu- 
lar beam epitaxy on n+ GaAs substrate. The undoped 
active layers consist of a 5.1nm quantum well sandwiched 
between 5nm Alo.33Gao.67 As barriers. Two 500nm GaAs 
layers, which are n-doped with Si to 10 18 cm~ 3 , pro- 
vide the top and bottom contacts to the double barrier 
structure. The barriers are separated from the contacts 
by undoped spacer layers to prevent Si segregation in 
the QD. Circular pillars of radii varying between 50nm 
and 1/xm were fabricated using electron-beam lithogra- 
phy and Ar/Ch reactive ion etching. A polyimide pla- 
narisation patterning was used for the bonding pads . 



Fig. |l|a shows the I-V characteristics of diodes as a 
function of the lateral size, at a temperature of AK . 
The voltage threshold shifts to higher bias when the ra- 
dius decreases and the current density decreases as well. 
This reflects the growing influence of the lateral depletion 
length. In the following we will focus only on the experi- 
mental data collected on the smallest pillar whose radius 
is lOOnm. A magnetic field is applied parallel to the cur- 
rent flow. Fig. [j]b shows the I-V curves at B — OT and 
9T at a temperature of 4K. A clear staircase-like struc- 
ture appears in the low voltage bias region, whose last 
steps are poorly resolved because of the sharp increase 
of the current. Similar staircase structures have been 
widely reported by different groups and originate cither 
from impurities or from QD's created by lateral confine- 
ment. Whenever a localized state is aligned with the 
emitter Fermi level, resonant tunneling is allowed and a 
plateau appears in the I-V characteristics. We argue in 
the following that, in our case, the plateaus are due to 
the discreteness of the N— electron states in the QD. 

Fig. |2| shows the dynamical conductance dl/dV of the 
pillar as a function of the voltage bias at three tempera- 
tures ranging from T — IK to T = 5K. From the tem- 
perature dependence of all the peaks labeled by |,0,H, X , 
and + in the figure we conclude that they correspond to 
the tunneling of electrons through QD states at the emit- 
ter Fermi level. In fact, taking into account the intrinsic 
width of the QD states, we fit the width at half maximum 
of the first conductance peak by using the formula: 

dl/dVcx / F(E - E F )—[T{E -aV)]dE (1) 

where F(E) is the Fermi distribution, Ep is the Fermi 
energy in the emitter, a is the energy-to- voltage con- 
version factor, and T(E) is a transmission probabil- 
ity with a Lorentzian shape. As a result we extract 
a ~ 0.31meV/mV, which is in good agreement with 
the thickness of the barriers and and the spacers. More- 
over, using this value of a, we check that the LO-phonon 
replica is at about 36 ± 5meV after the main resonance 
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peak (see fig. |l]b). Thus, a is fairly constant over the 
whole bias range. 

We performed magneto-tunneling measurements at 
T = 35mK applying a magnetic field along the direc- 
tion parallel to the current flow. Because the current has 
a sharp increase with the voltage close to the threshold, 
we cannot directly plot the conductance data for a quan- 
titative analysis. In fig. ||the derivative of the differential 
conductance, d 2 I/dV 2 is depicted in a gray-scale graph. 
This has no direct physical interpretation by itself, but it 
emphasizes the peaks of the conductance. The maximum 
of a conductance peak is represented in this diagram by 
the boundary line between a clear band and a dark band, 
in increasing the bias. The contour lines represent the 
dG/dV = OS.V" 1 values. 

This diagram also reveals fluctuations of the Local 
Density of State (LDoS) of the emitter Q. These fluc- 
tuations produce a fan-type pattern superimposed to the 
one due to the QD states, which becomes clearly visible 
at relatively high magnetic field. Indeed, for B > 2T, 
the LDoS fluctuations follow the behavior of the Landau 
levels in the emitter: the slope oc \dV\/dB of the fluc- 
tuations decreases with increasing B and corresponds to 
Landau levels which are deeply below the emitter Fermi 
level. We point out that their slope increases whenever a 
QD ground state is crossed. These slope changes are due 
to the progressive increase of the electron number in the 
QD. For B < 2T the Landau levels are destroyed by dis- 
order and the LDoS fluctuations are random. Although 
the LDoS fluctuations are clearly resolved, their current 
amplitude is weak compared to the one of the QD states 
and they produce very different features in the gray-scale 
plot. Thus, the two phenomena can be distinguished in 
almost the whole range of B and V. 

We focus now on the six clearly resolved boundary lines 
of the gray-scale, which correspond to the peaks labeled 
by an arrow in fig |[ We attribute them to the N-electron 
ground states with N=l,2,3,4,5 and 6. Let us assume a 
parabolic lateral confining potential V — ^m*ojor 2 and 
put the first peak at ~ 348mV in fig. ^ in correspondence 
with the single particle level E no of the Darwin-Fock 
spectrum H E nm = Hui(2n + \m\ + 1) + ^mhw c (with 
u> = \J lo 1 /4 + 0J 2 and ui c = eB/m*). Because we take 
a = 0.31meV/mV, we obtain Hloq = 9 ± lmeV. This 
corresponds to a radius r ~ ^j2h/ra*ujQ ~ 15 ± lnm. 
This is fairly larger than the expected radius of a single 
Si donor in GaAs, which is the most likely expected im- 
purity in the quantum well region. On the contrary this 
value of the radius is in very good agreement with the one 
extracted from the calculation of the lateral confining po- 
tential in the pillar. This calculation consists in solving 
the Poisson equation in the Thomas-Fermi approxima- 
tion for the electronic density in the contacts ||. The 
resulting band profile enables us to determine the energy 
levels of the QD by solving numerically the Schodinger 



equation. The single particle energy spectrum that is 
obtained is very well approximated by an harmonic po- 
tential with Hluq = 7 — 8meV, and the calculated radius 
of the ground state is equal to ~ 17nm. 

Despite the fact that the barriers of our structures have 
the same thickness, we have strong evidence that the sec- 
ond conductance peak (~ 368mV in fig. ||) originates 
from Coulomb blockade. The simplest way to include in- 
teraction is to assume that the two electron ground state 
is given by U% = 2E ofi + E C - We define as the n-th 
state of the N-electron spectrum, and Ec is the Coulomb 
interaction: 

E c = J\n^)t^£-^ d r^ (2) 

Since we measure the electro-chemical potential of the 
QD defined by fi(N) = — f/jy_i, we interpret the 
separation between the first and the second peak as Ec- 
There is a factor of two discrepancy, however, between 
the experimental value Eq = a.20mV ~ 6 — 7meV and 
the theoretical one estimated perturbatively from eq. || 
with single particle wavefunctions of the harmonic po- 
tential, E^ = e 2 y/m*u) /32iThe 2 ~ 14meV. This can be 
justified by taking into account the finite width of the 
QD as well as the screening of the electrons in the con- 
tacts. Following pi], we estimate the in plane separation 
of the electrons in the QD as lOnm, while the electrons 
in the contacts are at a distance A ~ 13nm (including 
the screening length in the contacts ~ 8nm). Therefore, 
the Coulomb integral should be cut off at separations 
r 2 — i"i > A, what reduces the bare theoretical Coulomb 
interaction of our device to the expected value. 

One more argument in favour of our interpretation 
is the observation that the separation between the two 
peaks increases slightly with B. Such an increase has 
already been observed jjj and is due to the squeezing of 
the electron orbits which causes an increase of Ec- 

At B = QT, the Darwin-Fock spectrum has degener- 
ate sets of states separated by Hujo, these shells being 
completely filled for N = 2,6, 12.... The degeneracy is 
completely lifted by the electron-electron interaction. We 
attribute the set of four lines in the voltage bias range 
420mV — 460mV to the N-electron ground states with 
N=3,4,5 and 6 which fill progressively the second shell. 
Keeping a constant in the whole bias range, the exper- 
imental separation between the two shells is ~ 15meV, 
which is about Ec + ft^oo, as expected in the Constant 
Interaction (CI) model. 

Fig. ^ is the result of an exact calculation of the en- 
ergy spectrum of the QD using a Lanczos algorithm, 
assuming a parabolic potential and a Coulomb interac- 
tion (huif) — 8meV,Ec = 6meV). The shell structure is 
clearly visible in fx(N), as well as the characteristic pair- 
ing of levels for particles with opposite spins. According 
to the CI model, the addition energy fi(N) — fi(N — 1) 
for each even electron should always be Ec- However 
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Eq is reduced to 5meV in the second shell, because of 
correlation effects |1 and we measure an even smaller 
value E^ p ~ 3 — 4meV in this shell. Again, we attribute 
the difference between E 1 ^ and E^ p mainly to the fact 
that the assumption of the QD as isolated in the Lanc- 
zos calculation is unphysical. It ignores the coupling of 
the dot with the electrodes. The occupation rate for the 
N > 4 electron ground states appears to be noticeably 
less than that of the first ground states, what implies the 
observed reduction of the Coulomb blockade. Indeed, in 
the real device, any increase of V lowers the collector 
barrier height and increases the energy window through 
which the tunneling occurs, hence enhancing transport 
through excited states. 

We conclude that there is a close correspondence be- 
tween the solid lines of fig. || and the six main lines of 
the experimentally observed conductance diagram and 
attribute these to tunneling while the dot is in its 
N— electron ground state, with increasing N. 

We claim that the additional temperature dependent 
features of the fig. 0, not discussed up to now, are due 
to the excited states of the N-electron spectra. Such 
separation of excited and ground states in a non-linear 
transport experiment implies that the collector barrier is 
slightly less transparent than the emitter one. This orig- 
inates probably from fluctuations of the nominal param- 
eters during the growth. This hypothesis is corroborated 
by the data in the negative bias, where the distinction 
between excited and ground states is not possible. These 
peaks appear to have a correspondence to analogous fea- 
tures in the calculated spectrum of fig. [|, lying between 
the first and the second shell. The 0, H and o-dotted lines 
are respectively the first, second and third one-electron 
excited states corresponding to the _Eo,-i,-Eo,i arL d Bo, -2 
Darwin-Fock states. The x line corresponds to the first 
excited 2 electrons state TJ\, with a total spin 5=1, 
referred to the single particle ground state energy 
(as all the U% in the following). Note that the singlet- 
triplet transition (t/^-f/jO occurs for B > 9T and is not 
observed. The + line describes the second 2-electrons ex- 
cited state Ul with S = 0. At B — 0T, C/f - rj 2 ° = fio> , 
so that we expect to observe the states in the following 
order when V increases: U®, U®, U\ (which is 2-fold de- 
generate), f/j, U$ and XJ\. 

In fig. U we have marked the experimental boundary 
lines that we attribute to these excited states. White 
symbols refer to lines that are clearly visible in our data. 
Black symbols indicate the expected position of these 
lines where the data are confused because of the LDoS 
fluctuations. Despite these uncertainties, the agreement 
between experimental and numerical results is satisfac- 
tory. We stress that, as in other experiments PJlQ], the 
degeneracy of J5q,±i is lifted at B = 0T. This is a con- 
sequence of the lack of inversion symmetry of the Zinc 
Blende structure. Moreover the current contribution of 
the .E?o,i washes out quickly with B as expected from the 



conservation of the angular momentum in 3D-0D tunnel- 
ing@Q. 

Let us finally comment on the origin of the main res- 
onance peak observed in the I-V curves (fig. 0b), like 
in large, non-laterally confined diodes involving tunnel- 
ing through 2D states. Such a peak appears only if the 
momentum perpendicular to the interfaces is conserved 
during the tunneling. This is not the case for transport 
through a QD where only the energy is conserved. The 
only way to explain the shape of the I-V characteristics 
is to add a continuum of states which arises from the 
excited states of the QD. 

In a simple non-interacting particles picture there are 
12 levels for tunneling (including the 4 ground states of 
the second shell) when the chemical potential is twice 
tkoQ. As the Coulomb interaction suppresses the degen- 
eracy, the mean level spacing is roughly ~ lmeV- which 
corresponds to our resolution in this experiment. The 
continuum becomes effective when /i ~ 3hu>o because of 
the 64 available states (including the 6 ground states of 
the third shell); as a consequence there is no way to ob- 
serve excited levels beyond the second shell of our QD 
structure. 

In conclusion, we have reported on non linear magneto- 
tunneling experiment conductance measurements carried 
out in a single QD. Our data are consistent with the spec- 
troscopy of few-electron states, even if LDoS fluctuations 
decrease the resolution of the experiment. We clearly 
observe a shell structure with a sizeable reduction of the 
Coulomb interaction in the second shell, probably due to 
screening. The first excited states for the one- and two- 
electron states can be identified. Finally we attribute the 
main I-V peak to the growing number of excited-states 
when the voltage bias is increased. Quantitative discrep- 
ancies with exact diagonalization of an isolated QD show 
that the role of the contacts in these non linear transport 
experiments cannot be ignored. 
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FIG. 1. (a) I-V curves at T = 4K for different pillars: 
curves refer to increasing radii (100,200,300, 350 and 1000 
nm), from bottom to top; (b) I-V curve of the smallest pillar 
(r = lOOnm). 
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FIG. 3. Gray-scale plots of the derivative of the differen- 
tial conductance = d 2 I/dV 2 , vs bias voltage and magnetic 
field (steps O.lmV and 0.02T) numerically obtained from the 
measured G(V) data, (white, dG/dV > 5/j.S/mV; black, 
dG/dV < —5/j.S/mV; dG/dV — corresponds to the con- 
tour lines. 
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FIG. 4. Chemical potential /i(iV) from the calculated en- 
ergy spectrum. The solid lines correspond to the dot in the 
ground states with N < 6; the dotted lines correspond to the 
dot in the lowest excited states for electron occupancy N = 1 
and 2 (see text for details). 



FIG. 2. Spectra of dl/dV vs V at T = 1, 2 and 5K. 
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